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Q\ . ABSTRACT 

Qh' We observed the anomalous X-ray pulsar 4U 0142+61 using the Proportional 

jy*l • Counter Array (PCA) aboard the Rossi X-ray Timing Explorer (RXTE) in 

■ March 1996. The pulse frequency was measured as v = 0.11510039(3) Hz 

with an upper limit of \v\ < 4 x 10~ 13 Hz s _1 upon the short term change 
>- ■ in frequency over the 4.6 day span of the observations. A compilation of 

all historical measurements showed an overall spin-down trend with slope 
z> = —3.0 ± 0.1 x 10~ 14 Hz s -1 . Searches for orbital modulations in pulse arrival 
times yielded an upper limit of a x sku ^ 0.26 lt-s (99% confidence) for the 
period range 70 s to 2.5 days. These limits combined with previous optical 
limits and evolutionary arguments suggest that 4U 0142+61 is probably not a 
~Q_f member of a binary system. 
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1. Introduction 

Most known accreting X-ray pulsars are members of high-mass X-ray binaries (HMXBs) 
containing an OB companion ( [van Paradijs 199~5| ; |Bildsten et al. 1997| ). Comparatively, 



the number of pulsars found in low-mass X-ray binaries (LMXBs) is quite small. Only 
five pulsars are known to be members of LMXBs: Her X-l, GX 1+4, GRO J1744-28, 4U 
1626-67, and SAX J1808.4-3658 = XTE J1808-369 QChakrabarty & Morgan TOB]) , which 



have known orbital periods or companions ( [Bildsten et al. 1997| and references therein.) In 
addition to the five known LMXB pulsars, there are 6 to 8 pulsars which are not members 
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of HMXBs: 4U 0142+61, IE 1048.1-5937, IE 2259+586, RX J1838.4-0301 (see |Mereghetti 
Belloni, fc Nasuti 1997] for an alternative interpretation of this system), IE 1841-045, 
1RXS J170849.0-400910 ( | Vasisht E Gotthelf 1997] and references therein), AX J1845-0258 
flGotthelf fc Vasisht 1998] ), and possibly RX J0720.4-3125 ( [Haberl et al. l997|) . These 
pulsars, called anomalous or braking X-ray pulsars ([van Paradijs, Taam, fc van den Heuvel 
|1995| ; [Mereghetti fc Stella 1995| ; |Ghosh, Angelini, fc White 1997| ), all have spin periods in 
the 5-11 s range, which is extremely narrow compared to the 69 ms-2.8 hr range of spin 
periods of HMXB pulsars. No optical counterparts of these systems have been detected 
to date, and 4U 0142+61 flSteinle et al. 1987| ; |White et aTT987| ; |Coe fc Pightling 1998|) , 
IE 1048.1-5937 ( |Mereghetti, Garaveo, fc Bignami 1991 ), IE 2259+586, RX J1838.4-0301 
( |Uoe fc Pightling 1998| and references therein), RX J0720.4-3125 flHaberl et al. 1997] ; 
[Motch fc Haberl 1998| ) have optical limits that rule out a high-mass companion. These 
pulsars are also characterized by the following properties: (1) a general spin down trend 
with time scales P/\P\ = 10 3 — 10 5 years, while HMXBs with similar periods show episodes 
of spin-up and spin-down ( [Bildsten et al. 1997|) ; (2) No evidence of orbital periodicity 
to date; (3) a low X-ray luminosity of 10 34 — 10 36 ergs s _1 that is quite constant on time 
scales of days to ~ 10 years; (4) a very soft X-ray spectrum usually well described by a 
combination of a blackbody with effective temperature ~ 0.3-0.4 keV and a photon power 
law (f(E) = AE-i) with photon index 7 = 3 - 4 QMereghetti, Stella, fc Israel 1998Q ; and 
(5) a relatively young age (^ 10 5 yr) as inferred from the low galactic scale height (~ 100 
pc r.m.s.) of the anomalous X-ray pulsars (AXPs) as a group ( |van Paradijs, Taam, fc van| 
den Heuvel 1995| ) and the apparent association of IE 2259+587 ( [Fahlman fc Gregory 198T| ), 
IE 1841-045 QVasisht fc Gotthelf 1997] ), and RX J1838.4-0301 ( ^chwentker T9"9l ) with 
supernova remnants. 



The persistent X-ray source 4U 0142+61 was discovered with Uhuru ([Forman et al 



1978| ). Pulsations at 8.7 s were discovered in 1984 August EXOSAT data ( [Israel, Mereghetti 



fc Stella 1994| ). A 25 minute modulation was also seen in these data ( |White et al. 1987|) , 
but was later found to be due to the nearby X-ray transient pulsar RX J1046. 9+6121 
( [Motch et al. 1991] ; [Hellier 1994| ). The ASCA 2-10 keV spectrum of 4U 0142+61 was 
well described by a blackbody with effective temperature 0.386 ± 0.005 keV plus a power 
law with photon index 3.67 ± 0.09 and an absorption column Nh = 9.5 ± 0.4 x 10 21 cm -2 
( |White et al. 1996| ). This spectrum was consistent with an earlier spectrum measured 
with EXOSAT when RX J0146.9+6121 was not present QWhite et al. 1987| ). The spin 
frequency of 4U0142+61 was measured with Einstein ( |White et al. 1996] ), EXOSAT ( [Israel 
[Mereghetti, fc Stella 1994] ), ROSAT ( Motch et al. 1991| ; [Hellier 1994| ), ASCA flWhite et 
al. 1996| ), and RXTE (this paper). The spin frequency history is shown in Figure |l|. All 
high significance historical measurements of the spin frequency follow the general spin-down 
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trend of z> = —3.1 ±0.1 x 10 14 Hz s l , which corresponds to a spin-down timescale z^/| z>| r 
123,000 yr. Three frequency measurements from 1985 EXOSAT post discovery data and 
1991 ROSAT data Qlsrael, Mereghetti, fc Stella 1994| ), which Israel et al. (1994) found to 



have high probabilities of chance occurrence, have been omitted from Figure |I| 

Israel et al. (1994) performed an orbital period search on ~ 12 hr of EXOSAT ME data 
from 1984 August 27-28. The search was carried out for periods from 430 s to 43,000 s, 
assuming a circular orbit. No significant orbital signature was found. Other EXOSAT 
measurements in 1985 November and December did not have sufficient time resolution to 
unambiguously detect 4U 0142+61. 

Here we present results from RXTE observations of 4U 0142+61, obtained in 1996 
March. The observation times and durations are listed in Table 1. An orbital period search 
is performed for periods between ~ 600 s and 4.6 days. Power spectra are presented and 
compared to white noise spectra with identical windowing. Upper limits are placed on the 
size of the allowed orbits for different trial periods. Next a modified algorithm is used to 
search for orbital periods shorter than ~ 600 s and to place upper limits on the size of 
allowed orbits. RXTE results are compared to previous results obtained with EXOSAT. 
Allowed companion types based upon these orbital limits and previous optical limits are 
discussed. 



2. Observations and Analyses 

The Proportional Counter Array (PCA) on RXTE consists of five xenon/methane 
multi-anode proportional counters sensitive to photons from 2-60 keV and has a total 
collecting area of 6500 cm 2 . The PCA is a collimated instrument with an approximately 
circular field-of-view with a FWHM of about 1° flJahoda et al. 1996Q . PCA observat ions of 



4U 0142+61 were performed on 1996 March 25 and 28-30. Details of these observations are 
given in Table 1. The field of view also contained the Be/X-ray pulsar RX J0146. 9+6121 
( |Motch et al. 1991| ; |Hellier 1994 [Haberl et al. 1998Q . For 1996 March 28-30, Good Xenon 



data (time and detector tagged events with 1 fis time resolution and full energy resolution) 
were used for our analysis. For 1996 March 25, Event Mode data, which have lower time 
resolution and similar energy resolution, were used. To maximize signal to noise ratio, 
data from only the top xenon layer of the PCA and from the energy range 3.7-9.2 keV 
were selected. The light curves were binned to 125 ms time resolution and the times were 
corrected to the solar system barycenter. 

The pulse frequency of 4U 0142+61 was determined to be z/ spin = 0.11510051(8) by an 
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epoch-folded search of data from 1996 March 25 and 28-30 for the frequency range 0.1139 
- 0.1165 Hz. Assuming the pulse frequency remained constant at this value across our 
observations, we divided the data into segments and determined phase offsets and intensities 
for each segment. A template profile was created by epoch- folding the first 80,000 seconds 
of data from March 28-30. The template profile is shown in Figure 0. The data were then 
divided into 250 second segments. In each segment the data were epoch-folded at the pulse 
frequency. The template and the pulse profiles from each segment were represented by a 
Fourier expansion in pulse phase. The profiles were limited to the first 3 Fourier coefficients, 
which contained 98.6% of the power in the template. These Fourier coefficients, a^, were 
calculated as 

1 N 

^ = ^T.R^ 3k/N (i) 

iV 3=1 

where k is the harmonic number, i = (— l) 1 / 2 , N = 16 is the number of phase bins in the 
pulse profile, and Rj is the count rate in the jth bin of the pulse profile. The pulse phase 
offset of each segment with respect to the template was calculated by cross-correlating each 
phase fold with the template. The cross-correlation is equivalent to a fit to 

3 In, — T P i2vk( l>t, | 2 

X 2 = E ' 2 ' (2) 

k=l a a k 

where is the Fourier coefficient of harmonic number k of the template profile, / is the 
relative intensity of the measured profile with respect to the template, and <fi is the phase 
offset of the measured profile relative to the template profile. The errors on the phase 
offsets were weighted according to the amount of data in each segment. Figure |3] shows the 
phase offsets with a best fit quadratic removed. A total of 7 points, which deviated from 
this fit by more than 0.25 cycles, were discarded. For 5 of these outliers, the 250 s bins 
were less than 1/3 full. The remaining two outliers, which both contained less than 200 
s of data and were offset by about 1/2 cycle, contained double-peaked pulse profiles. A 
quadratic fit to the remaining points yielded a best fit frequency of v = 0.11510039(3) Hz 
at MJD 50169.835 (consistent with our earlier epoch-folding results) and a 3 a upper limit 
of v £ 4 x 10~ 13 Hz s" 1 with x 2 /182 = 1-05. 

The phase offsets were searched for an orbital periodicity using a Lomb-Scargle 
periodogram ( press et al. 1992 and references therein) shown in Figure [|. The top panel is 



the Lomb-Scargle periodogram for the entire 4 day interval. The largest peak at « 586 s 
has a probability of chance coincidence of 65%. The center panel is the Lomb-Scargle 
periodogram for March 28-30 only. The largest peak at ~ 1680 s has a probability of chance 
coincidence of 78%. These probabilities were computed using Monte-Carlo simulations of 
10,000 white noise realizations. The phase offsets were simulated as white noise with zero 
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mean and standard deviations equal to the measurement errors on the phase offsets along 
with the actual time tags corresponding to the phase offsets. The cumulative probability 
distribution of the peak noise power was generated by retaining the power in the largest 
noise peak for each trial. The bottom panel is the Lomb-Scargle periodogram for white 
noise with mean and variance equal to that of the phase offsets with time tags corresponding 
to the 4 day interval. From the noise power, it is evident that much of the higher frequency 
variability in the power spectra is due to the windowing of the data. The Lomb-Scargle 
technique removes a constant from the data. Any variations in count rate that deviate from 
a constant result in low frequency variations in the power spectrum. This analysis revealed 
no significant orbital signatures. 

Previous upper limits, a x sini ^ 0.37 lt-s (430 s ^ P orb ^ 43000 s) estimated from 
EXOSAT data, by Israel et al. (1994), were calculated using the method of Van der Klis 
(1989) . This method assumes that the total power in a frequency bin is the sum of the 
signal power and the noise power. Vaughan et al. (1994) show that this assumption holds 
if a large number, n, of power spectra have been averaged together. However, if n is 
small or n — 1, as is the case for EXOSAT and RXTE measurements of 4U 0142+61, this 
assumption is incorrect. For small values of n, the total power is calculated from the vector 
sum of the Fourier amplitudes of the signal and noise. Vaughan et al. (1994) estimate 
that upper limits calculated using van der Klis (1989) will increase by at least 30% when 
calculated correctly. 

Upper limits for a x sini, the projected semi-major axis of the neutron star, were 
calculated under the assumption that signal and noise amplitudes are combined vectorially. 
No intrinsic spin frequency variations were detected, hence the pulse phase was given by 



were 0o is a constant phase and v is a constant frequency. The pulse emission time 
t em = t ssh — z where t ssb is the time corrected to the solar system bary center and z is the 
time delay due to binary motion. For a circular orbit the delay is given by, 



The phase offsets were modeled for an arbitrary epoch t (chosen as the mid-time of 
the data set) as 



= 00 + Vot' 



cm 



(3) 




(4) 



where P or b is the orbital period and T n / 2 is the epoch of 90° mean orbital longitude. 



0modcl = 00 + Az/ (t -t ) + A COS 
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for a fixed value of the trial orbital period P rb, where 0o was a constant offset in phase and 
Az/ is an offset in frequency. In this model, 

. . {A 2 + B 2 f 2 , R , 

a x smt = : — — _ (6) 

u smc (p— ) 



where z/ is the frequency used for epoch folding, sine x = and AT = 250 s, the bin size. 
The loss of sensitivity for short periods is accounted for by the sine (7rAT/P or b) windowing 
factor. The phase and frequency offsets, 0o and Av , were estimated for each orbital period, 
accounting for the loss of sensitivity at long orbital periods. For each trial orbital period, 
single parameter 99% confidence regions were constructed for A and B ( |Lampton, Margom 



fc Bowyer 1976| ). Since A = B = was not excluded for any trial orbital periods, an upper 



limit for a x sini was estimated from the maximum value of (A 2 + B 2 ) 1 ^ 2 , the point on 
the confidence region farthest from the origin. Upper limits on a x sin i at 99% confidence 
generated in this manner are shown in Figure |5|. Sensitivity is reduced for periods <^ 600 
seconds, for periods ^ 10 5 seconds, and at the RXTE orbital period (~ 96 min). An overall 
99% upper limit for a selected frequency range can be estimated from the largest upper 
limit within that range. For all periods in the range 587 s <; P or b <^ 2.5 days, a x sini ^ 0.27 
lt-s (99% confidence). If we exclude the region of reduced sensitivity around the RXTE 
orbital period, this limit reduces to 0.25 lt-s (99% confidence). 

To search for orbital periods shorter than about 600 seconds, a different method 
was employed. Well defined phase offsets and intensities could not be computed by 
cross-correlating profiles and the template for segments much shorter than 250 seconds. 
However, if the intensity was assumed to be known, phase offsets could be computed for 
much shorter intervals by approximating the cross-correlation function with a linear model. 
The intensity of 4U 0142+61 appeared to be reasonably constant on short time scales, so 
we assumed it was constant across 500 s intervals. Intervals 500 s long produced very good 
measurements of the phase offset and intensity, while allowing longer timescale intensity 
variations. Phase offsets and relative intensities were generated for each 500 second interval 
by the cross-correlation method described earlier. Next each 500 s interval was subdivided 
into shorter segments, each 2 pulse periods long (^ 17.4 seconds). Fourier coefficients were 
then fit to each 17.4 s segment by minimizing 

2 ^ {rj - [cq + SLi cos(27rfc0(t)) + v k sin(27r k<j){t))}} 2 

i=l U Ti 

where is count rate measurement i, N is the number of measurements in a 17.4 s segment, 
Co is a constant offset, k is the harmonic number, u k and Vk are the Fourier coefficients, 
and 4>{t) = 0o + vo(t — to) is the phase model. The Fourier coefficients then were fit by a 
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linearized model of the cross-correlation with x 2 given by 

2 >^ | a* - ^5oo(l + 27ri/cA(/>)t fc | 2 



x 2 = E 



fc=i I ""* 



2 



where = u fc — zt> fc , £ fe is the template, 7 50 o is the intensity for the corresponding 500 s 
interval, and A0 is the phase offset for each 17.4 s segment. The mean of the phase offsets 
for all 17.4 s segments within a 500 s interval corresponded to the phase offset computed for 
that 500 s interval using the cross-correlation method. The top panel of Figure |^ shows the 
Lomb-Scargle periodogram generated from the short interval phases. The largest peak at ~ 
87 s has a probability of chance coincidence of 18%. The bottom panel of Figure |6| shows 
the 99% confidence upper limits on a x sin % for fixed trial orbital periods. For the period 
range, 70 s ^ P or b ^ 610 s, a x sim <^ 0.26 lt-s at 99% confidence. 

EXOSAT ME data with 1 second time resolution from 1984 August were obtained from 
the High Energy Astrophysics Science Archive Research Center (HEASARC). Spacecraft 
position information was not easily obtainable, so the times were corrected only for 
motion of the Earth. The data in 1000 second segments were epoch-folded at the period 
-Pspin = 8.68723(4) s measured by Israel et al. (1994). A template profile was generated by 
epoch-folding the entire data set at the pulse period. The template and the pulse profiles 
from each segment were represented by a Fourier expansion in pulse phase. The pulse phase 
offset of each segment with respect to the template was calculated by cross-correlating each 
phase fold with the template. The phase offsets then were searched for an orbital periodicity 
using a Lomb-Scargle periodogram. No significant orbital signatures were found. Using the 
method described earlier, upper limits were placed on a x sini. For all periods in the range 
2050 s <z P or b <z 2 days, a x sini ^ 0.73 lt-s (99% confidence). This limit, generated under 
the assumption that signal and noise amplitudes are combined vectorially, is about twice 
the limit of a x sinz <^ 0.37 lt-s (99% confidence) obtained by Israel et al. (1994) . The latter 
was obtained under the assumption ([van der Klis 1989 ) that signal power and noise power 



may be added, which does not apply in this case. 



3. Discussion 

The limits on a x sini, in conjunction with an assumed inclination angle and an assumed 
neutron star mass of 1.4 Mq, can be recast to give a limit at each orbital period on the 
companion mass. These limits include no assumptions as to the nature of the companion. 
A plot of these limits, for inclination angles of 8, 30, and 90° versus orbital period is shown 
in Figure 0. 
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In a low-mass X-ray binary, the companion is expected to fill its Roche Lobe. Eggleton 
(1983) gave a useful expression for the Roche Lobe radius, given by 

R2 = 0-49g 2 / 3 

a 0.6g 2 /3 + i n (i + g i/3) W 

where a is the separation between the two stars, q = M 2 /M x is the mass ratio, M 2 is the 
companion mass, and M x is the neutron star mass. This expression can be combined with 
Kepler's third law to obtain the binary period as a function of the companion star's average 
density, p = 3M 2 /47ri?2 an d mass ratio q, given by 

/ - \ -1/2 / \ 1/2 
iW = 3 x 10 4 s ^ J (j^-j (0.6 + q~^ ln(l + fl*))*f* (10) 

where Penary is the binary period. If the structure of the Roche Lobe filling star is assumed, 
and thus a relation between its mass and radius, the binary period can be expressed as 
a unique function of the mass, M2, for an assumed value of M x = IAMq. Period-mass 
relations derived assuming a normal main sequence star, a helium main sequence star, a 
white dwarf flVerbunt fc van den Heuvel 1995Q and a period core mass relation for low 



mass 2Mq) giants flPhinney fc Kulkarni 19941) are listed in Table 2 and plotted in 



Figure [7|. The period core mass relation for low mass giants is valid for core masses of 
0.16Mq ^ M core ^ 0.45Mq, where the lower limit corresponds to the helium core mass at 
the end of main sequence evolution and the upper limit is the core mass at helium flash 
( IPhinney fc Kulkarni 1994| ). 

The maximum allowed mass of a Roche lobe filling hydrogen main sequence companion 
is ~ 0.25Mq for inclinations greater than 30° (99% confidence). This mass corresponds to 
an M5 or later type star. Larger masses are allowed if the inclination angle is quite small. 
For % = 8° the maximum mass is ~ 0.5Mq (99% confidence). The probability of observing 
a system with i < 8° by chance is 1%. If instead, a Roche lobe filling helium burning 
companion is assumed, higher masses (99% confidence) of ~ 0.65Mq for i > 30° and 2Mq 
for i = 8° are allowed. Low mass giants, similar to the companion to GRO J1744-28 ( [Finger 



et al. 19*96 ), with core masses of ~ 0.21Mq — 0.45Mq are allowed for i > 30° and core 



masses of ~ 0.18Mq — 0.45Mq are allowed for i > 8°. Period-mass relations predict white 
dwarf masses well below our upper limits for all inclinations. 



Optical (V £ 24, [Steinle et al. 1987| R £ 22.5, |White et al."T987|) and infra-red (J £ 



19.6 and K ^ 16.88, |Coe fc Pightling 1998|) observations failed to detect any counterpart 



within the ROSAT error circle ( [Hellier 1994Q . The optical limits can be used to further 
constrain the allowed companion masses based upon calculated optical magnitudes of a 
normal main sequence star, helium main sequence star, and the optical emission produced 
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by reprocessed X-rays in the accretion disc. The column density, Nh = 8 x 10 21 cm 2 
measured by White et al. (1996) , corresponds to an optical extinction of A v ~ 4.7 ( [Predehl 



fc Schmitt 1995| ), which implies a distance of ~ 5 kpc ( |Hakkila et al. 1997] ) . However, the 



dust scattering halo observed by ROSAT ( | White et al. 1 99(f ) was only half that predicted 



by Njj, suggesting that the absorbing material is clumped along the line of sight and that 
4U 0142+61 is most likely at a distance of < 5 kpc. A main sequence companion with a 
mass of ^ 0.8Mq (for any inclination angle) is compatible with the optical limits. 

The optical observations place much more stringent constraints on a giant companion 
with a helium core and on a helium main sequence companion. The minimum core mass 
for a giant star with a helium core is 0.16Mq, corresponding to the core mass when the 
star leaves main sequence. A giant of solar metalicity with a core mass of 0.16 Mq has 
a luminosity of L ~ Lq and an effective temperature of ~ 5000 K (|Phinney fc Kulkarni 



1994). At the upper limit distance of 5 kpc with absorption A v ~ 4.7 and a bolometric 



correction of -0.2 (|Zombeck 1990|) , we calculate an apparent magnitude of V ~ 23, for 



a giant star with a core mass of 0.16Mq. Hence a giant companion similar to that of 
GRO J1744-28 is not allowed by the optical limits. The minimum mass for a helium main 
sequence star is M ~ 0.3Mq (|Kippenhahn fc Weigert 1990|) . From plots in Kippenhahn & 
Weigert (1990), a 0.4Mq He star corresponds to a luminosity of ~ 5.6Lq and an effective 
temperature of ~ 30,000 K. At the upper limit distance of 5 kpc with absorption A v ~ 4.7 
and bolometric correction of -3 ( [Flower 1996| ), we calculate an apparent magnitude V ~ 24 



for a 0.4Mq helium star. Hence a helium star of a mass of <^ 0.4Mq is compatible with 
optical limits. 

In low-mass X-ray binaries, the optical emission is usually dominated by reprocessing 
of X-rays in the accretion disk. From fits to several low-mass X-ray binaries with known 
distances, van Paradijs and McClintock (1994) derived an empirical relationship between the 
absolute visual magnitude and the X-ray luminosity and orbital period. This relationship 
assumes that the companion is filling its Roche-lobe, the optical emission is dominated by 
the accretion disk, and that accretion disks are axially symmetric scaled up versions of a 
standard shaped toy model. No assumptions are made about the companion type. This 
relationship is given by 



M v = 1.57 - 2.27 log 



lhr/ \L 



>Edd 



where M v is the absolute visual luminosity, P is the binary period, L x is the X-ray 
luminosity, and L^dd = 2.5 x 10 38 erg s _1 is the Eddington luminosity. Using this relation 
along with estimates of the X-ray luminosity, L x = 7.2 x 10 34 d 2 pc flWhite et al. 199"^ ), 



absorption A v ~ 4.7 ([White et al. 1996|) , and apparent visual mag nitude V £ 24 QSteinle et| 



10 



al. 1987), we derived a relation between binary period P, and distance d. If the disk was 
fainter than the observed optical limit, then the allowed orbital periods are given by 



where 77 is an estimate of the scatter in the van Paradijs & McClintock relationship. From 
Figure 2 in Van Paradijs & McClintock (1994), we estimate a scatter in M v of ~ 0.75 
magnitudes in the data used to obtain equation (jTTf) . A decrease in M v of 1.5 magnitudes 
(twice the scatter), results in an increase in the maximum period by a factor of rj ~ 10. 
Thus, for the upper limit distance of 5 kpc (where L x = 2 x 10 36 erg s -1 , at which equation 
(|TT|) applies) this relation, including scatter, requires a period P <; 2300 s. Substituting this 
into the period- mass relationships in Table 2, results in allowable masses for only helium 
main sequence companions and white dwarf companions at 5 kpc. At distances <^ 3 kpc, 
which are more likely based upon ROSAT scattering measurements flWhite et al. 199~6|) , 
only white dwarf companions are allowed. 



We found no evidence (99% confidence) for orbital modulation in the pulse arrival 
times for orbital periods between 70 seconds and 2.5 days. Searches for orbital modulations 
in pulse arrival times yielded an upper limit of a x sku <^ 0.26 lt-s (99% confidence) for 
the period range 70 s to 2.5 days. Our limits on a x sini lead to 99% confidence dynamical 
limits for i ^ 30° of <^ 0.25Mq for normal main sequence companions, ^ 0.65Mq (99% 
confidence) for helium main sequence companions, and 0.21Mq ^ M corc ^ 0.45Mq (99% 
confidence) for <^ 2Mq giants with helium cores. Optical limits at 5 kpc allow masses of 
<^ 0.8Mq and <^ 0.4Mq for normal and helium main sequence companions respectively. 
Optical limits do not allow giant companions with helium cores. Optical and dynamical 
limits currently do not constrain white dwarf companions. 

The smooth spin-down (Figure |l|) observed in 4U 0142+61 is inconsistent with the 
random walk behavior expected for wind fed accreting pulsar (pildsten et al. 1997|) . 
Hence, 4U 0142+61 is unlikely to be a wind accretor. Long-term observations with BATSE 
( Bildsten et al. 1997| ) show that disk fed accreting pulsars switch between states of 
spin-up and spin-down with the magnitude of the torque in either state comparable to a 
characteristic torque. A pulsar subject to this torque will spin up (or down) at a rate of 




(12) 



4. 



Conclusion 



( IBildsten et afTT^ ) 




(13) 
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where M is the mass accretion rate and P sp i n is the pulsar spin period in seconds. The 
observed spin-down rate of i> = —3 x 10 -14 Hz s _1 corresponds to M ^ 9 x 10~ 12 Mq yr _1 
or a luminosity of L ^ 10 35 erg s" 1 which is consistent with the luminosity of 7.2d| pc x 10 34 
erg s _1 measured with ASCA ( |White et al. 199Tf ). Hence the observed long-term spin-down 
rate of 4U 0142+61 is consistent in behavior with the disk fed accreting pulsars. However, 
optical limits only allow accretion disks in systems with orbital periods <^ 2300 s. Orbital 
periods this small allow only helium main sequence and white dwarf companions at 
distances of 3-5 kpcs and only white dwarf companions at ^ 3 kpc,which are more likely 
based upon ROSAT scattering measurements ( [White et al. 1996| ). However, current 



evolutionary scenarios do not provide a mechanism for producing a Roche lobe filling 
white dwarf companion or a Roche lobe filling helium main sequence companion within the 
assumed ^ 10 5 yr age of the system. This age is based upon the low scale height of AXPs 
as a group ( [van Paradijs, Taam, fc van den Heuvel 1995| ) and ages inferred from supernova 
remnant associations of other AXPs. 

Optical, dynamical, and evolutionary limits argue against binarity in AXPs. Van 
Paradijs et al. (1995) proposed that AXPs are powered by accretion onto the neutron 
star from a circumstellar disc resulting from a common envelope evolution. It is not clear, 
however, why such an evolution would result in such a narrow range of pulse periods. The 
narrow period range is much more easily explained by the magnetar model (Thomson fc 



Duncan 1996| ). In this model, the neutron star surface is heated by the decay of the very 



strong magnetic field, producing X-rays. Heyl & Hernquist (1997) also propose a model 
consisting of a highly magnetized young neutron star surrounded by a thin envelope of 
hydrogen or helium. In this model, X-rays are produced by thermal emission. Pulsations 
are produced by a temperature gradient on the neutron star's surface induced by the strong 
magnetic field combined with limb darkening in the neutron star's atmosphere. Detailed 
long-term monitoring observations of the AXPs are needed to search for spin-up episodes 
which would suggest the system is an accretor, or a glitch accompanied by a soft-gamma 
repeater burst ( [l'homson fc Duncan 1996| ), which would suggest the magnetar model. 



This research made use of data obtained from the High Energy Astrophysics Science 
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Table 1: 


1996 RXTE Observations 


of 4U 0142+61 


Observation Start (TT) 


Observation Stop (TT) 


Observation Duration (ksec) 


March 25:12:31:04 


14:14:06 


5.533 


March 28:22:15:12 


March 29:00:02:06 


3.367 


March 29:00:11:28 


06:46:06 


13.987 


March 29:07:23:28 


11:34:06 


10.182 


March 29:22:18:40 


March 30:03:34:06 


11.104 



Table 2: Derived Mass-Orbital Period Relations for LMXBs a 

Companion Type Mass-Orbital Period Relation 

Giant with He core b P orb = 1.1 x 10 5 s (M corc /O.16M )/(g) c 

H Main Sequence P orb = 3 x 10 4 s (M c /M Q )f(q) 

He Main Sequence P orb = 3 x 10 3 s (M c /M & )f(q) 

White Dwarf P orb = 40 s (M /M c ) d 



"Adapated from Verbunt and van den Heuvel (1995) 

b Valid for O.16M £ M core £ O.45M flPhinncy fc Kulkarni 19941 ) 

c f(q) = (l + l/g)- 1/2 (0.6 + g- 2 / 3 ln(l + g 1 /3))3/2 

d For q < 0.8, (i.e. M c < 1.1M ), f{q) w 1 
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Fig. 1. — The historical pulse frequency history. The RXTE measurement is indicated by a 
filled circle. The dotted line is the best linear fit to the overall spin down v = — 3. 1±0. 1 x 1CT 14 
Hz s" 1 . 
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Fig. 2. — The mean-subtracted template pulse profile generated by epoch-folding the first 
80,000 s of 3.7-9.2 keV data from March 28-30 at a constant frequency, z/ spin =0. 11510041 Hz. 
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Fig. 3. — Phase offsets for 250 s intervals of data from 1996 March 25-30. A best fit 
quadratic has been subtracted. The time plotted is in days since 1996 March 25:07:12:00 
(MJD 50167.3). 
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Fig. 4. — Lomb-Scargle periodograms (normalized to trie data variance) tor the entire March 

25-30 data set (top panel), the March 28-30 dataset only (center panel), and white noise 

(bottom panel) with mean and variance equal to that of the phase offsets and time tags 

corresponding to the March 25-30 interval. The dotted lines indicate the 99% confidence 

level for a detection. 
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Fig. 5. — Upper limits to a x sini (99% confidence) for a circular orbit at a trial orbital period. 
Sensitivity is reduced for both the longest and shortest trial orbital periods. Removal of data 
containing Earth occultations also reduces sensitivity at the RXTE orbital period. Curves 
of constant mass function are shown. The dotted line denotes our overall 99% confidence 
upper limit of a x sini <^ 0.26 lt-s. 
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Fig. 7. — Companion mass limits calculated using the 99% confidence limits on a x sim from 
Figures [5] & || assuming a 1.4 Mq neutron star and 8, 30, 90° inclination angles. Also shown 
are the period-mass relations for a normal main sequence star, helium burning star, and white 
dwarf and a period core mass relation for a giant star. 



